home *** CD-ROM | disk | FTP | other *** search
/ Pascal Super Library / Pascal Super Library (CW International)(1997).bin / MATH / NRPAS13 / FTEST.DEM < prev    next >
Text File  |  1991-04-29  |  1KB  |  58 lines
  1. PROGRAM d13r7(input,output);
  2. (* driver for routine FTEST *)
  3. CONST
  4.    npts=1000;   (* npts must be > mpts for FTEST *)
  5.    mpts=500;
  6.    eps=0.04;
  7.    nval=10;
  8. TYPE
  9.    gln1array = ARRAY [1..npts] OF real;
  10.    gln2array = ARRAY [1..mpts] OF real;
  11.    glnparray = gln1array;
  12. VAR
  13.    glinext,glinextp : integer;
  14.    glma : ARRAY [1..55] OF real;
  15.    gliset : integer;
  16.    glgset : real;
  17.    f,factor,prob,vrnce : real;
  18.    i,idum,j : integer;
  19.    data1 : gln1array;
  20.    data2 : gln2array;
  21.  
  22. (*$I MODFILE.PAS *)
  23. (*$I RAN3.PAS *)
  24.  
  25. (*$I GASDEV.PAS *)
  26.  
  27. (*$I GAMMLN.PAS *)
  28.  
  29. (*$I BETACF.PAS *)
  30.  
  31. (*$I BETAI.PAS *)
  32.  
  33. (*$I AVEVAR.PAS *)
  34.  
  35. (*$I FTEST.PAS *)
  36.  
  37. BEGIN
  38. (* generate two gaussian distributions with 
  39. different variances *)
  40.    gliset := 0;
  41.    idum := -144;
  42.    writeln;
  43.    writeln('Variance 1 = ':16,1.0:5:2);
  44.    writeln('Variance 2':13,'Ratio':11,'Probability':16);
  45.    FOR i := 1 to (nval+1) DO BEGIN
  46.       FOR j := 1 to npts DO BEGIN
  47.          data1[j] := gasdev(idum)
  48.       END;
  49.       vrnce := 1.0+(i-1)*eps;
  50.       factor := sqrt(vrnce);
  51.       FOR j := 1 to mpts DO BEGIN
  52.          data2[j] := factor*gasdev(idum)
  53.       END;
  54.       ftest(data1,npts,data2,mpts,f,prob);
  55.       writeln(vrnce:11:4,f:13:4,prob:13:4)
  56.    END
  57. END.
  58.